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ABSTRACT 

We consider a class of geometric facility location problems in which 
the goal is to determine a set X of disks given by their centers (tj) 
and radii (ry) that cover a given set of demand points K C R^ at 
the smallest possible cost. We consider cost functions of the form 
Y,if{rj), where /(r) = r" is the cost of transmission to radius r. 
Special cases arise for a = 1 (sum of radii) and a = 2 (total area); 
power consumption models in wireless network design often use an 
exponent a > 2. Different scenarios arise according to possible re- 
strictions on the transmission centers tj, which may be constrained 
to belong to a given discrete set or to lie on a line, etc. 

We obtain several new results, including (a) exact and approx- 
imation algorithms for selecting transmission points tj on a given 
line in order to cover demand points y C R^; (b) approximation 
algorithms (and an algebraic intractability result) for selecting an 
optimal line on which to place transmission points to cover Y\ (c) 
a proof of NP-hardness for a discrete set of transmission points in 
R- and any fixed a > 1 ; and (d) a polynomial-time approximation 
scheme for the problem of computing a minimum cost covering 
tour (MCCT), in which the total cost is a linear combination of the 
transmission cost for the set of disks and the length of a tour/path 
that connects the centers of the disks. 
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1. INTRODUCTION 

The problem. We study a geometric optimization problem that 
arises in wireless network design, as well as in robotics and vari- 
ous facility location problems. The task is to select a number of 
locations tj for the base station antennas {servers), and assign a 
transmission range rj to each tj, in order that each p, e Y for a 
given set y = {pi, .. .,p„} of « demand points {clients) is covered. 
We say that client p; is covered if and only if p,- is within range of 
some transmission point tj., i.e., d{tj.,pi) ^ rj.. The resulting cost 
per server is some known function /, such as f{r) = r". The goal is 
to minimize the total cost, L;/(0')' "^^'^ ^1' placements of at most 
k servers that cover the set Y of clients. 

In the context of modeling the energy required for wireless trans- 
mission, it is common to assume a superlinear (a > 1) dependence 
of the cost on the radius; in fact, physically accurate simulation 
often requires superquadratic dependence (a > 2). A quadratic de- 
pendence {a = 2) models the total area of the served region, an ob- 
jective arising in some applications. A linear dependence (a = 1) 
is sometimes assumed, as in Lev-Tov and Peleg [19], who study 
the base station coverage problem, minimizing the sum of radii. 
The linear case is important to study not only in order to simplify 
the problem and gain insight into the general problem, but also 
to address those settings in which the linear cost model naturally 
arises [10,21]. For example, the model may be appropriate for a 
system with a narrow-angle beam whose direction can either ro- 
tate continuously or adapt to the needs of the network. Another 
motivation for us comes from robotics, in which a robot is to map 
or scan an environment with a laser scanner [13, 14]: For a fixed 
spatial resolution of the desired map, the time it takes to scan a cir- 
cle coffesponds to the number of points on the perimeter, i.e., is 
proportional to the radius. 



Our problem is a type of clustering problem, recently named 
min-size k-clustering by Bilo et al. [7]. Clustering problems tend to 
be NP-hard, so most efforts, including ours, are aimed at devising 
an approximation algorithm or a polynomial-time approximation 
scheme (PTAS). 

We also introduce a new problem, which we call minimum cost 
covering tour (MCCT), in which we combine the problem of find- 
ing a short tour and placing covering disks centered along it. The 
objective is to minimize a linear combination of the tour length 
and the transmission/covering costs. The problem arises in the au- 
tonomous robot scanning problem [13,14], where the covering cost 
is linear in the radii of the disks, and the overall objective is to 
minimize the total time of acquisition (a linear combination of dis- 
tance travelled and sum of scan radii). Another motivation is the 
distribution of a valuable or sensitive resource: There is a trade- 
off between the cost of broadcasting from a central location (thus 
wasting transmission or risking interception) and the cost of trav- 
elling to broadcast more locally, thereby reducing broadcast costs 
but incurring travel costs. 

Location Constraints. In the absence of constraints on the server 
locations, it may be optimal to place one server at each demand 
point. Thus, we generally set an upper bound, k, on the number of 
servers, or we restrict the possible locations of the servers. Here, 
we consider two cases of location constraints: 

(1) Servers are restricted to lie in a discrete set {/i , . . . ,fm}; or 

(2) Servers are constrained to lie on a line (which may be fuUy 
specified, or may be selected by the optimization). 

Our results. We provide a number of new results, some improv- 
ing previous work, some giving the first results of their kind. 

In the discrete case studied by Lev-Tov and Peleg [19], and Bilo 
et al. [7], we give improved results. For the discrete ID problem 
where F C K, we improve their 4-approximation to a linear-time 
3-approximation by using a "Closest Center with Growth" (CCG) 



algorithm, and, as an alternative to the previous 0{(n - 



algo- 



rithm [19], we give a near-linear- time 2-approximation that uses a 
"Greedy Growth" (GG) algorithm. Unfortunately, we cannot ex- 
tend our proofs to the l^D problem. Intuitively, greedy growth 
works as follows: start with a disk with center at each server, each 
disk of radius zero; among all clients, find one that requires the 
least radial disk growth to capture it; repeat until all clients are cov- 
ered. Note that for a ^ 2 the 2D variants of the problem are already 
proved to be NP-Hard and to have a PTAS [7]. 

In the general 2D case with clients K C K , we strengthen the 
hardness result of Bilo et al. [7] by showing that the discrete prob- 
lem is already hard for any superlinear cost function, i.e., f{r) = r" 
with a > 1. Furthermore, we generalize the min-size clustering 
problem in two new directions. On the one hand, we consider less 
restrictive server placement policies. For instance, if we only re- 
strict the servers to lie on a given fixed line, we give a dynamic 
programming algorithm that solves the problem exactly, in time 
0{n log n) for any Lp metric in the linear cost case, and in time 
0{n'^ log«) in the case of superlinear non-decreasing cost functions. 
For simple approximations, our algorithm "Square Greedy" (SG) 
gives in time 0{nlogn) a 3-approximation to the square covering 
problem with any linear or superlinear cost function. A small varia- 
tion, "Square Greedy with Growth" (SGG), gives a 2-approximation 
for a linear cost function, also in time 0(;ilog«). The results are 
also valid for covering by Lp disks for any p, but with correspond- 
ingly coarser approximation factors. 

A practical example in which servers are restricted to lie along a 
line is that of a highway that cuts through a piece of land, and the 



server locations are restricted to lie along the highway. The line lo- 
cation problem arises when one not only needs to locate the servers, 
but also needs to select an optimal corridor for the placement of the 
highway. Other relevant examples may include devices powered by 
a microwave or laser beam lining up along the beam. 

If the servers are restricted to lie on a horizontal line, but the 
location of this line may be chosen freely, then we show that the 
exact optimal position (with a = 1) is not computable by radicals, 
using an approach similar to that of Bajaj [5, 6] in addressing the 
unsolvability of the Fermat- Weber problem. On the positive side, 
we give a fully polynomial-time approximation scheme (FPTAS) 
requiring time 0((n-'/e)log;i) if a = 1 and time 0{{n^ /£) log n) if 
a> 1. 

For servers on an unrestricted line, of any slope, and a = I, we 
give 0(l)-approximations (4-approximation in 0(;!^log«) time, or 
8\/2-approximation in 0{n log n) time) and an FPTAS requiring 
time 0{{n^ /e^) log n). 

We give the first algorithmic results for the new problem, mini- 
mum cost covering tour (MCCT), which we introduce. Given a set 
y C R- of « clients, our goal is to determine a polygonal tour T and 
a set X of k disks of radii rj centered on T that cover Y while min- 
imizing the cost length(r) +C^r". Our results are for a = 1. The 
ratio C represents the relative cost of touring versus transmitting. 
We show that MCCT is NP-hard if C is part of the input. At one 
extreme, if C is small then the optimum solution is a single server 
placed at the circumcenter of Y (we can show this to be the case 
for C ^ 4). At the other extreme (if C very large), the optimum 
solution is a TSP among the clients. For any fixed value of C > 4, 
we present a PTAS for MCCT, based on a novel extension of the 
m-guillotine methods of [20]. 

Related work. There is a vast family of clustering problems, 
among which are the k-center problem in which one minimizes 
maxj rj, the k-median problem in which one minimizes Y,id{pi,tji), 
and the k-clustering problem in which one minimizes the maximum 
over all clusters of the sum of pairwise distances between points in 
that cluster. For the geometric instances of these related clustering 
problems, refer to the survey by Agarwal and Sharir [I]. When k is 
fixed, the optimal solution can be found in time 0{n ) using brute 
force. In the plane, one of the only results for the min-size cluster- 
ing problem is a small improvement for fc = 2 by Hershberger [17], 
in subquadratic time 0(n /loglogn). Approximation algorithms 
and schemes have been proposed, particularly for geometric in- 
stances of these problems (e.g., [4]). Clustering for minimizing the 
sum of radii was studied for points in metric spaces by Charikar 
and Panigrahy [9], who present an 0(l)-approximation algorithm 
using at most k clusters. 

For the linear-cost model (a = 1), our problem has been consid- 
ered recently by Lev-Tov and Peleg [19] who give an 0((/7 + m) ) 
algorithm when the clients and servers all lie on a given line (the 
ID problem), and a linear-time 4-approximation in that case. They 
also give a PTAS for the two-dimensional case when the clients and 
servers can lie anywhere in the plane. Bilo et al. [7] show that the 
problem is NP-hard in the plane for the case f(r) = r", a ^ 2, either 
when the sets X and Y are given and k is left unspecified {k = n), 
or when k is fixed but then X = Y . They give a PTAS for the linear 
cost case (a = 1) and a slightly more involved PTAS for a more 
general problem in which the cost function is superlinear, there are 
fixed additive costs associated with each transmission server and 
there is a bound k on the number of servers. 

There are many problems dealing with covering a set of clients 
by disks of given radius. Hochbaum and Maass [18] give a PTAS 
for covering with a minimum number of disks of fixed radius, where 



the disk centers can be taken anywhere in the plane. They intro- 
duce a "grid- shifting technique," which is used and extended by 
Erlebach et al. [12]. Lev-Tov and Peleg [19] and Bilo et al. [7] 
extend the method further in obtaining their PTAS results for the 
discrete version of our problem. 

When a discrete set X of potential server locations is given, Gon- 
zalez [16] addresses the problem of maximizing the number of 
covered clients while minimizing the number of servers supply- 
ing them, and he gives a PTAS for such problems with constraints 
such as bounded distance between any two chosen servers. In [8], 
a polynomial-time constant approximation is obtained for choosing 
a subset of minimum size that covers a set of points among a set of 
candidate disks (the radii can be different but the candidate disks 
must be given). 

The closest work to our combined tour/transmission cost (MCCT) 
is the work on covering tours: the "lawn mower" problem [2], and 
the TSP with neighborhoods [3,11], each of which has been shown 
to be NP-hard and has been solved with various approximation al- 
gorithms. In contrast to the MCCT we study, the radius of the 
"mower" or the radius of the neighborhoods to be visited is speci- 
fied in advance. 

2. SCENARIO (1): SERVER LOCATIONS RE- 
STRICTED TO A DISCRETE SET 

2.1 The one-dimensional discrete problem with 
linear cost 

Consider the case of m fixed server locations X = {ti,...,!,,,}, n 
client locations Y = {pi , ...,p„}, and a linear (a = I) cost function, 
with clients and servers all located along a fixed line. Without loss 
of generality, we may assume that X and Y are sorted in the same 
direction, at an extra cost of 0{{n + m) log{n + m)). Lev-Tov and 
Peleg [19] give an 0{{n + m)^) dynamic programming algorithm 
for finding an exact solution. Bilo et al. [7] show that the problem is 
solvable in polynomial time for any value of a by reducing it to an 
integer linear program with a totally unimodular constraint matrix. 
The complexities of these algorithms, while polynomial, is high. 
Lev-Tov and Peleg also give a simple "closest center" algorithm 
(CC) that gives a linear-time 4-approximation. We improve to a 
3-approximation in linear time, and a 2-approximation in 0(m + 
nlogm) time. 

We now describe an algorithm which also runs in linear time, but 
achieves an approximation factor of 3. 

Closest Center with Growth (CCG) Algorithm: Process the clients 
{p\,...,Pn} from left to right keeping track of the rightmost ex- 
tending disk. Let cos denote the rightmost point of the rightmost 
extending disk, and let R denote the radius of this disk. (In fact the 
rightmost extending disk will always be the last disk placed.) If cos 
is equal to, or to the right of the next client processed, p,, then pi 
is already covered so ignore it and proceed to the next client. If p; 
is not yet covered, consider the distance of p, to cor compared with 
the distance of p; to its closest center f,. If the distance of p/ to 
co« is less than or equal to the distance of p, to its closest center fj, 
then grow the rightmost extending disk just enough to capture p;. 
Otherwise use the disk centered at /)■ of radius |p; — //| to cover p;. 

Lemma 1. For a = 1, CCG yields a 3-approximation to OPT 
in 0{n-\-m) time. 

The proof is similar to that of the next lemma, and omitted in this 
version. 



If we consider a single disk D with clients pi and pR on the 
left and right edges of £), associated centers xi, xr at distances 
respectively radius(D)— e to the left and radius(D)— e to the right, 
along with a dense set of clients in the left hand half of D we see 
that 3 is the best possible constant for CCG. 

Finally we offer an algorithm that achieves a 2-approximation 
but runs in time 0(m + « log m). 

Greedy Growth (GG) Algorithm: Start with a disk with center 
at each server all of radius zero. Now, amongst all clients, find 
the one which requires the least radial disk growth to capture it. 
Repeat until all clients are covered. An efficient implementation 
uses a priority queue to determine the client that should be captured 
next. One can set up the priority queue in 0{m) time. Note that the 
priority queue will never have more than 2m elements, and that 
each Pi eventually gets captured, either from the right or from the 
left. Each capture can be done in time O(logra) for a total running 
time of 0(ra + ;ilogm). 

Lemma 2. For a = 1, GG yields a 2-approximation to OPT in 
0{m-\-n\ogm) time. 

Proof. Define intervals /, as follows: when capturing a client p, 
from a server tj whose current radius (prior to capture) is ry, let 
Jj = {tj + rj,pi] if Pi > tj, and 7/ = [pi,tj — rj) otherwise. Our first 
trivial yet crucial observation is that 7; n Jj^ = & if i ^ k. Also note 
that the sum of the lengths of the 7, is equal to the sum of the radii 
in the GG cover 

Consider now a fixed disk D in OPT, centered at f/j, and the list 
of intervals y, whose p; is inside D. As before, at most one such 7, 
extends outward to the right from the right edge of D. If so, call it 
Jr, and define Ji symmetrically. If 7^ exists, it cannot extend more 
than radius(D) to the right of D. Let A, = length (7s). We argue that 
there is an interval of length A. in D, to the right of to, which is free 
of 7,'s. It follows that there is at most radius(D) worth of segments 
to the right of f/). Of course, this is also true if 7r does not exist. By 
symmetry, there is also radius(£)) worth of segments to the left of 
tD, whether 7/, exists or not, yielding the claimed 2-approximation. 

Assume 7r exists. Then the algorithm successively extends Jg 
by growth to the left up to some maximum point (possibly stop- 
ping right at ps). Since the growth could have been induced by 
clients to the right of 7r, that maximum point is not necessarily 
a client. There is, however, some client inside D that is captured 
last in this process. This client p, (possibly pr) cannot be within 
A. of t£), since otherwise it would have been captured prior to the 
construction of 7r. 

If there is no client between f/j and p, we are done, since then 
there could be no interval 7,t in between. Thus consider the client 
p,_i just to the left of p;. Suppose d(p,_i,p;) ^ X. Then, if p,_i 
is eventually captured from the left, we would have the region be- 
tween p,_i and Pi free of 7jr.'s and be done. On the other hand, if 
Pi-i is captured from the right, it must be captured by a server be- 
tween p,_ 1 and Pi, and that server is at least A. to the left of p; since 
otherwise p; would be captured by that server prior to ps. This 
leaves the distance from the server to p, free of 7i's. 

Hence the only case of concern is if c?(p,_i,p,) < X. Clearly 
p,_i must not have been captured at the time when ps is captured 
since otherwise p, would have been captured before ps, contradict- 
ing the assumption that p, is captured by growth leftward from ps. 
Similarly, there cannot be a server between p,_i and p,, since oth- 
erwise both p,_i and p,- would be captured before pR. Together 
with the definition of p,, this implies that p,_i is captured from the 
left. Therefore, to the left of p,_ i , there must be one or more inter- 
vals {7;,} whose length is at least A that are constructed before p,_i 



is captured. Similarly, to the right of p/, there must be some one or 
more intervals {Jr } whose length is at least A., constructed before 
Pi is captured. However, either the last Jj. is placed before the last 
Ji- or vice versa. In the first case, there are no X length obstructions 
left in the left-hand subproblem, so pj_ i will be covered, and with 
X length obstructions remaining in the right subproblem, p,- will be 
captured by growth rightward. The second case is symmetrical to 
the first. In either case we have a contradiction. □ 

To see that the factor 2 is tight, just consider servers at — 2 + e, 
and 2 — e and clients at — 1 and 1 . 

2.2 Hardness of the two-dimensional discrete 
problem with superlinear cost 

In 2D, we sketch an NP-hardness proof, for any a > 1 . This strength- 
ens the NP-hardness proof of [7], which only works in the case 
a ^2. 

Theorem 3. For any a fixed a > 1, let the cost function of 
a circle of radius r be f{r) = r". Then it is NP-hard to decide 
whether a discrete set of n clients in the plane, and a discrete set 
ofm potential transmission points allow a cheap set of circles that 
covers all demand points. 

Proof (sketch). Let / be an instance of PLANAR 3SAT, and let G/ 
be the corresponding variable-clause incidence graph. After choos- 
ing a suitable layout of this planar graph, resulting in integer vari- 
ables with coordinates bounded by a polynomial in the size of G/ 
for all vertices and edges, we replace each the vertex represent- 
ing any particular variable by a closed loop, using the basic idea 
shown in the left of FigureQ this allows two fundamentally differ- 
ent ways of covering those points cheaply (using the "odd" or the 
"even" circles), representing the two truth assignments. For each 
edge from a vertex to a variable, we attach a similar chain of points 
that connects the variable loop to the clause gadget; the parity of 
covering a variable loop necessarily assigns a parity to all incident 
chains. Note that choosing sufficiently fine chains guarantees that 
no large circles can be used, as the overall weight of all circles in 
a cheap solution will be less than 1. (It is straightforward to see 
that for any fixed a > 1 , this can be achieved by choosing coordi- 
nates that are polynomial in the size of G/, with the exponent being 
0(l/(a-l)).) 



client points 
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Figure 1. (Left) The switch structure of a variable gadget. Note how there 
are two fundamentally different ways to cover all points cheaply (Right) The 
structure of a clause gadget. One small circle is needed for picking up the 
client point at the center of the gadget. 

For the clauses choose a hexagonal arrangement as shown in the 
right of Figure Q There is one central point that must be covered 
somehow; again, a > 1 guarantees that it is cheaper to do this from 
a nearby transmission point, rather than increasing the size of a 
circle belonging to a chain gadget. 



Now it is straightforward to see that there is a cheap cover, using 
only the forced circles, iff the truth assignment corresponding to 
the covering of variabe loops assures that each clause has at least 
one satisfying variable. O 

3. SCENARIO (2): SERVER LOCATIONS 
RESTRICTED TO A LINE 

3.1 Servers along a fixed horizontal line 

3.1.1 Exact solutions 

Suppose that the servers are required to lie on a fixed horizontal 
line, which we take without loss of generality to be the x-axis. Such 
a restriction could arise naturally (e.g., the servers must be con- 
nected to a power line, must lie on a highway, or in the main corri- 
dor in a building). In addition, this case must be solved first before 
attempting to solve the more general problem — along a polygonal 
curve. 

In this section, we describe dynamic programming algorithms to 
compute a set of server points of minimum total cost. For notational 
convenience, we assume that the clients Y are indexed in left-to- 
right order Without loss of generality, we also assume that all the 
clients lie on or above the x-axis, and that no two clients have the 
same x-coordinate. (If a client p, lies directly above another client 
Pj, then any circle enclosing p; also encloses pj, so we can remove 
Pj from Y without changing the optimal cover) 

Let us call a circle C pinned if it is the leftmost smallest axis- 
centered circle enclosing some fixed subset of clients. Equivalently, 
a circle is pinned if it is the leftmost smallest circle passing through 
a chosen client or a chosen pair of clients. Under any Lp metric, 
there are at most 0{n ) pinned circles. As long as the cost func- 
tion / is non-decreasing, there is a minimum-cost cover consisting 
entirely of pinned circles. 

Linear Cost. If the cost function / is linear (or sublinear), we 
easily observe that the circles in any optimum solution must have 
disjoint interiors. (If two axis-centered circles of radius r; and rj 
intersect, they lie in a larger axis-centered circle of radius at most 
r; + rj.) In this case, we can give a straightforward dynamic pro- 
gramming algorithm that computes the optimum solution under any 
Lp metric. 

The algorithm given in Figure |2| (left) finds the minimum-cost 
cover by disjoint pinned circles, where distance is measured using 
any Lp metric. We call the rightmost point enclosed by any pinned 
circle C the owner of C. 

If we use brute force to compute the extreme points enclosed by 
each pinned circle and to test whether any points lie directly above 
a pinned circle, this algorithm runs in 0{n ) time. With some more 
work, however, we can improve the running time by nearly a linear 
factor 

This improvement is easiest in the Loo metric, in which circles 
are axis-aligned squares. Each point p,- is the owner of exactly / 
pinned squares: the unique axis-centered square with p, in the up- 
per right corner, and for each point pj to the left of p,, the leftmost 
smallest axis-centered square with p, and pj on its boundary. We 
can easily compute all these squares, as well as the leftmost point 
enclosed by each one, in 0(/log i) time. (To simplify the algorithm, 
we can actually ignore any pinned square whose owner does not lie 
on its right edge.) If we preprocess P into a priority search tree in 
0{nlogn) time, we can test in 0{logn) time whether any client lies 
directly above a horizontal line. The overall running time is now 
0(n~logn). 



For any other Lp metric, we can compute the extreme points en- 
closed by all 0{n ) pinned circles in 0{n ) time using the follow- 
ing duality transformation. If C is a circle centered at (x,0) with 
radius r, let C* be the point {x,r). For each client pi, let p* = 
{C* \C is centered on the .\--axis and p; e C}, and let Y* = {p* \ p, e 
Y}. We easily verify that each set p* is an infinite x-monotone 
curve. (Specifically, in the Euclidean metric, the dual curves are 
hyperbolas with asymptotes of slope ±1.) Moreover, any two dual 
curves p^ and p*: intersect exactly once; i.e., Y* is a set of pseudo- 
lines. Thus, we can compute the arrangement of Y* in 0{n~) time. 
For each pinned circle C, the dual point C* is either one of the 
clients p, or a vertex of the arrangement of dual curves Y* . A cir- 
cle C encloses a client pi if and only if the dual point C* lies on 
or above the dual curve p*. After we compute the dual arrange- 
ment, it is straightforward to compute the leftmost and rightmost 
dual curves below every vertex in 0{n ) time by depth-first search. 

Finally, to test efficiently whether any points lie directly above 
an axis-centered (Lp) circle, we can use the following two-level 
data structure. The first level is a binary search tree over the x- 
coordinates of Y. Each internal node v in this tree corresponds to a 
canonical vertical slab 5,, containing a subset p,, of the clients. For 
each node v, we partition the j-axis into intervals by intersecting 
it with the furthest-point Voronoi diagram of py, in 0(|p,.| log|p,.|) 
time. To test whether any points lie above a circle, we first find a 
set of 0(log;i) disjoint canonical slabs that exactly cover the circle, 
and then for each slab 5,, in this set, we find the furthest neighbor 
in p,. of the center of the circle by binary search. The region above 
the circle is empty if and only if all 0(log«) furthest neighbors 
are inside the circle. Finally, we can reduce the overall cost of the 
query from 0(log n) to O(logM) using fractional cascading. The 
total preprocessing time is 0{nlog n). 

Theorem 4. Given n clients in the plane, we can compute in 
0(n log«) time a covering by circles (in any fixed Lp metric) cen- 
tered on the X-axis, such that the sum of the radii is minimized. 

Superlinear Cost. A similar dynamic programming algorithm 
computes the optimal covering under any superlinear (in fact, any 
non-decreasing) cost function /. As in the previous section, our al- 
gorithm works for any Lp metric. For the moment, we will assume 
that p is finite. 

Although two circles in the optimal cover need not be disjoint, 
they cannot overlap too much. Clearly, no two circles in the opti- 
mal cover are nested, since the smaller circle would be redundant. 
Moreover, the highest point (or apex) of any circle in the optimal 
cover must lie outside all the other circles. If one circle A contains 
the apex of a smaller circle B, then the lune 6\ A is completely con- 
tained in an even smaller circle C whose apex is the highest point 
in the lune; it follows that A and B cannot both be in the optimal 
cover. See Figure|3la). 




Figure 3. (a) The apex of each circle in the optimal cover lies outside the 
other circles, (b) The points Yc He in the shaded region, (c) If A and C are 
adjacent circles in the optimal covering, the shaded region B(A,C) is empty. 

To compute the optimal cover of y, it suffices to consider sub- 
problems of the following form. For each pinned circle C, let Y^ 



denote the set of clients outside C and to the left of its center; see 
Figure|3jb). Then for each pinned circle C, we have cost{Yc) = 
minj^{ f(radius{A)) -^cost{YA)), where the minimum is taken over 
all pinned circles A satisfying the following conditions: (1) The 
center of A is left of the center of C; (2) the apex of A is outside 
C; (3) the apex of C is outside A; and (4) A encloses every point in 
Yc\Ya. The last condition is equivalent to there being no clients 
inside the region B{A,C) bounded by the x-axis, the circles A and 
C, and vertical lines through the apices of A and C; see Figure|3jc). 

Our dynamic programming algorithm (Figure |2| (right)) consid- 
ers the pinned circles C] ,C2, . . . ,Cp in left to right order by their 
centers; that is, the center of C/ is left of the center of Cj whenever 
i < j. To simplify notation, let Yj = Y^.. For convenience, we add 
two circles Co and Cp+i of radius zero, centered far to the left and 
right of Y, respectively, so that Yq = and Yp^i = Y. 

Implementing everything using brute force, we obtain a run- 
ning time of 0{n^). However, we can improve the running time 
to 0{n log n) using the two-level data structure described in the 
previous section, together with a priority search tree. The region 
B(Ci,Cj) can be partitioned into two or three three-sided regions, 
each bounded by two vertical lines and either a circular arc or the x- 
axis. We can test each three-sided region for emptiness in O(logn) 
time. 

Theorem 5. Let f : E+ -^ M. be a fixed non-decreasing cost 
fimction. Given n clients in the plane, we can compute in 0{n logn) 
time a covering by circles (in any fixed Lp metric) centered on the 
X-axis, such that the sum of the costs of the circles is minimized. 

The algorithm is essentially unchanged in the Loo metric, except 
now we define the apex of a square to be its upper right comer. 
It is easy to show that there is an optimal square cover in which 
no square contains the apex of any other square. Equivalently, we 
can assume without loss of generality that if two squares in the 
optimal cover overlap, the larger square is on the left. To compute 
the optimal cover, it suffices to consider subsets Yc of points either 
directly above or to the right of each pinned square C. For any 
two squares A and C, the region B(A,C) is now either a three-sided 
rectangle or the union of two three-sided rectangles, so we can use 
a simple priority search tree instead of our two-level data structure 
to test whether B{A,C) is empty in 0(log;7) time. 

However, one further observation does improve the running time 
by a linear factor: Without loss of generality, the rightmost box in 
the optimal cover of Y(j has the rightmost point of Yq on its right 
edge. Thus, there are at most n candidate boxes Cj to test in the 
inner loop; we can easily enumerate these candidates in 0{n) time. 

Theorem 6. Let f : R+ ^ M foe a fixed non-decreasing cost 
function. Given n clients in the plane, we can compute in 0{n logn) 
time a covering by axis-aligned squares centered on the x-axis, 
such that the sum of the costs of the squares is minimized. 

3.1.2 Fast and simple solutions 

In this section we describe simple and inexpensive algorithms that 
achieve constant factor approximations for finding a minimum-cost 
cover with disks centered along a fixed horizontal line L, using any 
Lp metric. The main idea for the proofs of this section is to asso- 
ciate with a given disk D in OPT, a set of disks in the approximate 
solution and argue that the set of associated disks cannot be more 
than a given constant factor cover of Z), in terms of cumulative edge 
length, cumulative area, and so forth. 



MinSumOfRadiusCircleCover(I') : 
for every pinned circle C 

find tlie leftmost and rightmost points enclosed by C 
Co.rt[0]^0 
for (' ^ 1 to n 

Cost [('] ^- oo 

for each pinned circle C owned by pi 

if no points in P lie directly above C 
Pj <— leftmost point enclosed by C 
Cost[i] <— min{Co.v?[(], Cost[j — 1] + radius{C))} 
return Corf fn] 



MlNSUPERLINEARCOSTClRCLECOVER(r,/) : 


sort the pinned circles from left to right by their centers 


Cost[0] ^0 




for _/ ^ 1 to p + 1 




Cost [j] ^- CO 




for (■ ^ 1 to 7 — 1 




if C, and Cj exclude each other's apices | 


and B(C,,C;) 


is empty 


Costlj] 


- mm{CostlJ], Cost[{\+ f{mdius(Ci)))} 


return Cost[p+ 1] 





Figure 2. The dynamic programming algorithm: Left: linear cost; Right: superlinear cost function. 



As in the previous section, the case of Loo metric is the eas- 
iest to handle. By equivalence of all the Lp metrics, constant- 
factor c-approximations for squares will extend to constant-factor 
c'-approximations for Lp disks. 

Square Greedy Cover Algorithm (SG): Process the client points 
in order of decreasing distance from the line L. Find the farthest 
point p\ from L\ cover p\ with a square Si exactly of the same 
height as pi centered at the projection of p\ on L. Remove all 
points covered by 5] from further consideration and recurse, find- 
ing the next farthest point from L and so forth. In the case where 
two points are precisely the same distance from L, break ties arbi- 
trarily. 

Obviously, SG computes a valid covering of y by construction. We 
begin the analysis with a simple observation. 

Lemma 7. In the SG covering, any point in the plane {not nec- 
essarily a client) cannot be covered by more than two boxes. 

Proof. Suppose 5; and Sj are two squares placed during the run- 
ning of SG and that / < j so that 5; was placed before Sj. Then 
Si cannot contain the center point of Sj since then Sj would not 
have had the opportunity to be placed, and similarly Sj cannot con- 
tain the center point of Si. Now consider a point p e 5; fl Sj. If p 
were covered by a third square S^ then either one of {S,,5y} would 
contain the center of S^, or Sk would contain the center of one of 
{S/jSy}, neither of which is possible. □ 

Theorem 8. Given a set Y of n clients in the plane and any 
Ot ^ 1, SG computes in time 0(«logn) a covering of Y by axis- 
aligned squares centered on the x-axis whose cost is at most three 
times the optimal. 

Proof. Let Y = {pi,.. ..p,,} and consider a square S in OPT. 
We consider those squares {5, } selected by SG corresponding to 
points {pi : Pi e 5}, see Figure |4] and argue that these squares 




Figure 4. Squares of the SG algorithm inside a square of the optimal solu- 
tion. 



cannot have more than three times the total edge length of 5. The 



same will then follow for all of SG and all of OPT. The argument, 
without modification, covers the case of cost measured in terms 
of the sum of edge length raised to an arbitrary positive exponent 
a> 1. 

Arguing as in Lemma0it is easy to see that at most two boxes 
Si associated with points p; e 5 processed by SG actually protrude 
outside of 5, one on the left and one on the right. Denote by r the 
total horizontal length of these protruding parts of squares, then 
r ^ .?, the side length of S, since the side length of each protruding 
square is at most s and at most half of each square is protruding. 

Because of Lemma |7| the total horizontal length of all nonpro- 
truding parts of the squares 5; is at most 2s, consequently all points 
covered by 5 in OPT are covered by a set of squares 5, in SG whose 
total (horizontal) edge length Y.j^i is at most 3^. 

For exponents a > 1 observe that Y,j ^i < 35 and < .?, < .? for 
all i imphes that Eyi",- " < 3.?". 

To analyze the running time of the algorithm we need some more 
details about the data structures used: Initially, sort the points by 
x-coordinate and separately by distance from the line L in time 
0(n\ogn) and process the points in order of decreasing distance 
from L. As the point p,- at distance rf; from L is processed, we 
throw away points which are within horizontal distance d, from 
Pi. This takes time 0(log« + fei) time where A;,- is the number of 
points within di from p,-. Since we do this up to n times with 
k\ H h fc/ =n the total running time is 0(/7 log n). □ 

For the linear cost function, it is easy to modify the SG algorithm 
to get a 2-approximation algorithm. 

Square Greedy with Growth Algorithm (SGG): Process the points 
as in SG. However, if capturing a point p; by a square 5; would re- 
sult in an overlap with already existing square Sj then, rather than 
placing 5/, grow Sj just enough to capture p/, keeping the vertical 
edge furthest from pi at the same point on L. If placing 5/ would 
overlap two squares, grow the one which requires the smallest edge 
extension. Break ties arbitrarily. 

A proof somewhat similar to that of Lemma|2|shows that: 

Theorem 9. Given n clients in the plane, SGG computes in 
time 0{n\ogn) a covering by axis-aligned squares centered on the 
X-axis whose cumulative edge length is at most twice the optimal. 

Proof. As we process points pi using SGG, attribute to each point 
Pi a line segment .v/ along L as follows. If processing p,- resulted 
in the placement of a square 5; centered at the projection of pi in 
L then attribute to p, the projection on L of a horizontal edge of 
5/ (Case 1). If, on the other hand, processing of p; resulted in the 
growing of a prior square Sj to just capture p,-, attribute to p/ the 
projection on L of the portion of the horizontal edge of the ex- 
panded 5/ needed to capture p,- (Case 2). (This amount is at most 
the distance of p,- to L since otherwise p; would have been fallen 



into case 1.) We must show that the lengths of the segments is no 
more than twice the edge lengths of squares in OPT. 

It suffices to show that for any square 5 in OPT, the segments s, 
associated with points p, e 5 processed by SGG cannot have total 
edge length which exceeds twice the edge length s of 5. 

To see this observe that the sum of the lengths of those Sj lying 
completely inside S does not exceed .? since they are nonoverlap- 
ping. In addition, each of the parts of the at most two segments 
protruding from 5 can have length at most s/2, in case 1 for the 
same reason as in the SG algorithm, in case 2 since the total length 
of the segment is at most s/2. 

In order to make SGG efficient, we proceed as in SG. In addi- 
tion, we maintain a balanced binary search tree containing the x- 
coordinates of the vertical sides of the squares already constructed. 
For each new point pi to be processed we locate its x-coordinate 
within this structure to obtain its neighboring squares and to decide 
whether case 1 or case 2 applies. This can be done in time 0(log;i) 
just as adding a new square in case 1 or updating an existing square 
in case 2. Removing points covered by the new or updated square 
is done as in SG, so that the total runtime remains 0{nlogn). □ 

Unlike SG, SGG is not a constant factor approximation for area. 
Consider n consecutive points at height 1 separated one from the 
next by distance of 1 + e. Processing the points left to right using 
SGG covers all points with one square of edge length n + {n~ l)z, 
and so area 0{n ), while covering all points with n overlapping 
squares each of edge length 2, uses total area 4«. 

Finally, extending these results from squares to disks in any Lp 
metric is not difficult. Enclosing each square in the algorithm by 
an Lp disk leads to an approximation factor 3c^ for GG and 2c^ 
for SGG, where c = p"'''. In particular, for L2 disks, this yields a 
2\/2-approximation for a = I and a 4-approximation for a = 2. 

3.2 Finding the best axis-parallel line 

When the horizontal line £ is not given but its orientation is fixed, 
we first prove that finding the best line, even for a = 1, is uncom- 
putable, then in this linear case give a simple approximation, and 
finally a PTAS. 

3.2.1 A hardness result - uncomputability by radicals 

Our approach is similar to the approach used by Bajaj on the un- 
solvability of the Fermat- Weber problem and other geometric opti- 
mization problems [5, 6]. 

Theorem 10. Let c{t) = ^.tr, denote the minimum cost of a 
cover whose centers lie on the line of equation y = t. There exists a 
set Y of clients such that, iftQ is the value that minimizes c{t), then 
f() is uncomputable by radicals. 

The proof proceeds by exhibiting such a point set and showing 
by differentiating c{t) that tQ is the root of a polynomial which is 
proven not to be solvable by radicals. 

The following definitions and facts can be found in a standard 
abstract algebra reference; see, for example, Rotman [22]. A poly- 
nomial with rational coefficients is solvable by radicals if its roots 
can be expressed using rational numbers, the field operations, and 
taking kxh roots. The splitting field of a polynomial f{x) over the 
field of rationals Q is the smallest subfield of the complex numbers 
containing all of the roots of /(x). The Galois group of a poly- 
nomial f{x) with respect to the coefficient field Q is the group of 
automorphisms of the splitting field that leave Q fixed. If the Ga- 
lois group of /(x) over Q is a symmetric group on five or more 
elements, then f{x) is not solvable by radicals over Q. 



Consider the following set of points:{(3,4), (—3,-2), (102,2), 
(98,— 2), (200,— 2)}. By exhaustive case analysis, we can show 
that the optimal solution must consist of one circle through the first 
two points, a second circle through the next two points, and a third 
circle touching the last point, and the optimal horizontal line must 
lie in the range —2 <y<2. For a given value of v in this range, the 
cost of the best cover is 



c{y): 



2(,y- 1)2 + 18 + \/2y2 + 8 + (2->' 



Therefore, in order to find the best horizontal line, we must mini- 
mize c(>'). Setting the derivative to zero, we obtain the equation 



c'{y) 



2{y-i) 



2v 



V2(y-1)2 + 18 V2>'2 + 8 



1=0. 



We easily verify that c"(y) is always positive. The minimum value 
c{y) ^ 8.3327196 is attained at y ^ 1.4024709, which is a root of 
the following polynomial: 

/(y) = 1024 + 5 12y - 1600y^ + 1536y^ - 960y'' 
+368y5 - 172y'' + 28y^ - 7y'^ . 

Using the computational system GAP [15], we compute that the 
Galois group of /(y) is the symmetric group Sg, so the polynomial 
is not solvable by radicals. 

3.2.2 Fast and simple constant-factor approximations 

The simple constant factor approximations for a fixed line can be 
extended to the case of approximations to the optimal solution on 
an arbitrary axis-parallel line with the same constant factors, though 
with a multiplicative factor of 0{n^) increase in running time. 

3.2.3 An FPTAS for finding the best horizontal line 

We begin with the case a = 1. Let d denote the distance between 
the highest and lowest point. Clearly, d/2 ^ OPT $; nd/2. Partition 
the horizontal strip of height d that covers the points into 2n/e hor- 
izontal strips, each of height 5 = d^/2n, using 2n/z — 1 regularly- 
spaced horizontal lines, £/. For each line (i, we run the exact dy- 
namic programming algorithm, and keep the best among these so- 
lutions. Consider the line, I* , that contains OPT. We can shift line 
l* to the nearest i-,, while increasing the radius of each disk of OPT 
by at most 5, to obtain a covering of the points by disks centered 
on some ii; the total increase in cost is at most 5n = rfe/2 ^ e 
OPT. Thus, our algorithm computes a (1 + e) -approximation in 
time 0((;i-'/e)log«). 

In order to generalize this result to the case a > 1, let us write 
PSEUDO-OPT for the lowest cost of a solution on any of the hor- 
izontal lines £,, SHIFT for the result of shifting OPT to the closest 
of these lines, and r\ , ..., r,„ for the radii of the optimal set of disks. 
For an arbitrary power a > 1, we have 

m 

PSEUDO-OPT < SHIFT < ^ (r, + 5)" 

1=1 

m m 

< X-f + 8ai:(r,+8)«-i 

< OPT(l + 5a22'='-'«/d). 

The last line uses 8<d,ri<d and OPT > (rf/2)". Choosing 5 = 
Zd/{a2?-'^^^n) gives the desired (1 +e)-approximation. 

Together with the results from previous sections we have: 



Theorem 11. Given n clients in the plane and a fixed a ^ 1, 
there exists an FPTAS for finding an optimally positioned horizon- 
tal line and a minimum-cost covering by disks centered on that line. 
It runs in time 0{{n /e) logn) for the linear cost case (a = I) and 
0{{n^ /z)\ogn) for a> 1. 

3.3 Approximating the best line - 
any orientation 

Finally, we sketch approximation results for selecting the best line 
whose orientation is not given. We give both a constant factor ap- 
proximation and a PTAS for the linear cost case (a = 1). 

3. 3. 1 Fast and simple constant-factor approximations 

Given a line £, we say that a set D of disks D\,... ,£)j. is (.-centered 
if the centers of every disk Cj in D belongs to £. Recall that the 
cost of © is the sum of all its radii. 

Lemma 12. Given k ^ i, a line £, an l-centered set 1) of k 
disks that cover Y, and any point pQ on £, there exist p' eY and an 
£' -centered set 'D'ofk disks that cover Y, where (! is the line that 
joins pq and p', such that the cost ofl)' is at most 2" times the cost 
ofD. 

Proof. We will assume without loss of generality that £ is the x- 
axis, po is the origin and that no other point in Y lies on the y-axis. 
The latter restriction can easily be enforced by a small perturbation. 
Let the coordinates of p; be x; and v;, and let m,- denote the slope 
yi/xi of the line £i for 1 ^ i ^ n. First, we reorder Y so that \mi\ ^ 
. . . ^ |m„|. In what follows we assume that xi > and y\ ^ 0. The 
other cases can be treated analogously. 

For each disk Dj = D{tj,rj) in 2) , we construct a disk /)'■ whose 
radius is / = 2rj and center f'- is obtained from tj by rotating it 
around the origin counterclockwise by an angle tan^'(rai). The 
set 1)' of k disks thus defined is ^'-centered, where / = {(x,y) e 
R^ I y = raix} and pi e i". To see that d' covers Y, simply observe 
that d{tj,t':) ^ rj for all 1 ^ j ^ fc and apply the triangle inequality: 
any point in Dj must be at distance at most 2rj of f'-. The cost of 
this new solution is clearly at most 2" times that of O in the linear 
cost case. □ 

By a double appplication of this lemma, first about an arbitrary 
Po yielding a point p' = p,-, then about p; yielding another p' = pj, 
it is immediate that any i'-centered cover of Y can be transformed 
into an C;.j -centered cover whose cost is increased at most four- 
fold, where £;.,■ is the line joining p,- and py. By computing (exactly 
or approximately) the optimal set of disks for all 0{n ) lines de- 
fined by two different points of Y, we conclude: 

Theorem 13. Given n clients in the plane and a fixed a ^ 1, 
in 0{n logn) time, we can find a collinear set of disks that cover P 
at cost at most A'^OPT, and for a = 1, in 0{iv' logn) time, we can 
find a collinear set of disks that cover P at cost at most %\f2.0PT . 

3.3.2 A PTAS for finding the best line with 
unconstrained orientation 

We now prove that finding the best line with unconstrained orienta- 
tion and a minimum-cost covering with disks whose centers are on 
that line admits a PTAS. 

Theorem 14. Let Y he a set ofn clients in the plane that can 
be covered by an optimal collinear set of disks at linear cost OPT 
{i.e., a = I), and e > 0. In 0{{n'^ /E^)\ogn) time, we can find a 
collinear set of disks that cover Y at cost at most ( 1 + e) OPT. 



Proof. Let // be a strip of minimal width h that contains Y . Using 
a rotating calipers approach, H can be computed in 0{n\ogn) time. 
If h = 0, we can conclude that OPT = and we are done. 

Otherwise, we can assume wlog that H is horizontal and that its 
center line is the x-axis. Let R denote the smallest enclosing axis- 
parallel rectangle i? of y, w its width, and h its height. Then h <w 
and, moreover, h/2 ^ OPT. Let £* be the optimal line. 

We now distinguish two cases: 

Case 1. w ^ 2h: Observe that both vertical sides vi,V2 of R 
contain a point of Y. Therefore, (* must have distance at most 
OPT to vi and V2. A straightforward calculation shows that then £* 
must intersect the lines £i and £2 extending I'l and V2 at a distance 
of at most 40PT from the x-axis. 

The idea is now to put points on those parts of £1 and £2 which 
are equally spaced at distance 6 = zOPT/n. Then we consider all 
lines passing through one of these points on £1 and one on £2. For 
each such line we find the optimal covering of P by circles centered 
on it using the algorithm of Theorem |4| and give out the best one 
as an approximation for the optimum. 

Observe, that there is one of the lines checked, (, whose inter- 
section points with (i and £2 are at distance at most 8/2 from the 
ones of £*. Elementary geometric considerations show that to any 
point p in £* closest to some point of P there is a point in £ within 
distance at most 5. Consequently, to any circle of radius r of the 
optimal covering centered on £*, there is a circle on I of radius r + 5 
covering the same set of points (or more). Thus, / has a covering 
that differs by at most n8 from the optimal one. By the choice of 5 
we have a 1 + e-approximation to the optimum. 

Observe, that we chose 0{0PT/8) = 0{n/z) points on £1 and 
£2, so we are checking 0{{n/e) ) lines. For each of them, we apply 
the algorithm of Theorem|4|which has runtime O (n^ log n) yielding 
a total runtime of 0{n /e logn). 

Case 2. w < 2h: In this case the optimal line £* can have a 
steeper slope and even be vertical. Of course, it must intersect R 
and we expand R to a cocentric rectangle R' such that the foot- 
point of any point in Y on £* must lie inside R'. An easy geomet- 
ric consideration shows that extending the width of R by h and its 
height by w will suffice, so R' is a square of side length w + h. 
Then we put equally spaced points of distance 5 = zOPT/n on the 
whole boundary of ^', apply the algorithm of Theorem|4]to all lines 
passing through any two of these points, and return the one giving 
the smallest covering as an approximation to the optimum. The 
same consideration as in the first case shows that this is indeed a 
(1 + e)-approximation. Since the length of the boundary of R' is 
4{w + h) ^ 12/7 $; 240PT, we obtain the desired runtime in this 
case, as well. 

For both cases it remains to show how to obtain a suitable value 
of 5, since we do not know the value of OPT. Since any value 
below OPT suffices, we simply run a constant factor c approxima- 
tion algorithm of Theorem^and take 1 /c times the value it returns 
instead of OPT in the definition of 5. D 

4. MINIMUM-COST COVERING TOURS 

We now consider the minimum cost covering tour (MCCT) prob- 
lem: Given k'^ 1 and a set y = {pi , . . . ,p,i} of n clients, determine 
a cover of Y by (at most) k disks centered at X = {fi , . . . ,/j.} with 
radii ry and a tour T visiting A", such that the cost length(r) +CX1 ''" 
is minimized. We refer to the tour T, together with the disks cen- 
tered on X, as a covering tour of Y. Our results are for the case of 
linear transmission costs (a = 1). We first show a weak hardness 
result, then characterize the solution for C ^ 4, and finally give a 
PTAS for a fixed C> 4. 



4.1 A hardness result 



4.3 The case C> 4: A PTAS 



We prove the NP-hardness of MCCT where C is also part of the 
input. Note that this does not prove the NP-hardness of MCCT 
where C is a fixed constant, which is the problem for which we give 
a PTAS below. Note also that C appears in the run time exponent 
of that PTAS, and so the PTAS no longer runs in polynomial time 
if C is not a fixed constant. 

Theorem 15. MCCT with linear cost is NP-hard if the ratio C 
is part of the input. 

Proof (sketch). We show a reduction from HAMILTON CYCLE 
IN GRID GRAPHS. Given a set of n points on a grid, we construct 
an instance of MCCT in which each of the given points is a client. 
We set C to be larger than 2n. We claim that the grid graph has a 
Hamilton cycle if and only if there is a tour T visiting a set of disk 
centers with radii r,- whose cost is at most n. 

Clearly a Hamilton cycle in the grid graph yields a tour of cost 
n with each client contained in a disk of radius centered at that 
point. 

Conversely, suppose we have a tour whose cost (length plus sum 
of radii) is at most n. Note that no two clients can be contained 
in a single disk, as such a disk must have radius at least 0.5, and 
thus its contribution to the cost C • r, > 2« • 0.5 = n contrary to our 
assumption. Next we want to show that each disk in an optimal 
solution is centered at the client it covers. Suppose this is not the 
case, there is some client j which is covered by a disk centered 
at Cj ^ j. Let the distance between client j and the center of the 
disk covering it be d. Now consider an alternate feasible solution 
in which the tour visits Cj then j then back to cy, covering / with a 
disk of radius 0. No other client is affected by this change, as the 
disk only covers point j. The cost of the new solution is the cost of 
the original (optimal) solution +2d — Cd as we add 2d to the length 
of the tour, but decrease C^r, by Cd. Since C > 2 the new solution 
is better than the original optimal solution, a contradition. □ 

4.2 The case C ^ 4: The exact solution is 
a single circle 

Theorem 16. In the plane, with a cost function of length{T) + 
CY,fi cind C ^ 4, the minimum-cost solution is to broadcast to all 
clients from the circumcenter of the client locations and no tour 
cost. 

The proof rests on the following elementary geometry lemma 
(whose proof is omitted here). 

Lemma 17. For three points p, q and r in the plane, such that 
the triangle pqr contains its own circumcenter, the length of a trip 
from p to q to r and back to p is at least 4r where r is the circum- 
radius of the points. 

Proof of Theorem^] Let r{X) and r{Y) denote the minimum 
radius of a circle enclosing X orY, respectively. Let T be a covering 
tour of K, X C r be the set of disk centers and rj their radii. Finally, 
letrmax =maxjrj. 

By the triangle inequality, Lemma FTTl implies that the length(r) Js 
4r{X). Since the tour visits all the centers in X and the disks cen- 
tered at X cover Y, we have r{Y) ^ r{X) + rmax- By definition, the 
cost of T is length(r) +C£ ,■ rj, which by the observation above is 
at least 4r{X)+CY.jrj > 4r(X) -t-Crmax- The assumption C < 4 
then implies that it l3e at least C{r{X) + rnjax) ^ Cr{Y), which is 
the cost of covering by a single disk with a zero-length tour. □ 



We distinguish between two cases for the choice of transmission 
points: they may either be arbitrary points in the plane (selected by 
the algorithm) or they may be constrained to lie within a discrete 
set T of candidate locations. 

The constant C specifies the relative weight associated with the 
two parts of the cost function - the length of the tour, and the sum 
of the disk radii. If C is very small (C ^ 4), then the solution is to 
cover the set Y using a single disk (the minimum enclosing disk), 
and a corresponding tour of length (the singleton point that is 
the center of the disk). If C is very large, then the priority is to 
minimize the sum of the radii of the k disks. Thus, the solution 
is to compute a covering of y by fc disks that minimizes the sum 
of radii (as in [19]), and then link the resulting disk centers with 
a traveling salesman tour (TSP). (In the case that k ^ n, the disks 
in the covering will be of radius 0, and the problem becomes that 
of computing a TSP tour on Y.) Note that our algorithm gives an 
alternative to the Lev-Tov and Peleg PTAS [19] for coverage alone. 

Our algorithm is based on applying the ni-guillotine method [20], 
appropriately adapted to take into account the cost function and 
coverage constraint.' We need several definitions; we largely fol- 
low the notation of [20]. Let G = (V.E) be an embedding of a 
connected planar graph, of total Euclidean edge-length L. Let O be 
a set of disks centered at each vertex v of G of radius r,,. We refer 
to the pair (G, 2) ) as a covering network if the union Uvgvflr of the 
disks covers the clients Y. We can assume without loss of general- 
ity that G is restricted to the unit square B, i.e., Uggge C int{B). 

Our algorithm relies on there being a polynomial-size set of can- 
didate locations for the transmission points that will serve as the 
vertices of the covering tour we compute. In the case that a set T 
of candidate points is given, this is no issue; however, in the case 
that the transmission points are arbitrary, we appeal to the follow- 
ing grid-rounding lemma (proved in the full paper). 

Lemma 18. One can perturb any covering network {G,1)) to 
have its vertices all at grid points on a regular grid of spacing 
5 = 0{£ ■ diam{S)/n), while increasing the total cost by at most a 
factor of [l +e). 

An axis-aligned rectangle, W C i?, is called a window; rectan- 
gle W will correspond to a subproblem in a dynamic programming 
algorithm. An axis-parallel line I that intersects W is called a cut. 

For a covering network with edge set E and a set of disks © , we 
say that (£,©) satisfies the m-guillotine property with respect to 
window W if either (I) all clients K C W lie within disks of 2) that 
intersect the boundary of W; or (2) there exists a cut £ with certain 
properties (an m-good cut with respect to W) that splits W into W\ 
and Wi, and [E.V) recursively satisfies the ra-guillotine property 
with respect to both Wi and W2. Due to the lack of space, we cannot 
give the full definition of an ra-good cut (see the full paper). 

The crux of the method is a structural theorem, which shows how 
to convert any covering network (G, D ) into another covering net- 
work (G', d'), such that the new graph G' satisfies the nj-guillotine 
property, and that the total cost of the new instance (G',d') is 
at most 0{{L + CR)/m) times greater than the original instance 
(G, ©), where L is the total edge length of G and R the sum of 
the radii of ©. The construction is recursive: at each stage, we 
show that there exists a cut with respect to the current window W 
(which initially is the unit square B), such that we can "afford" (by 
means of a charging scheme) to add short horizontal/vertical edges 

'The "m" in this section refers to a parameter, which is 0(1/e), not the 
number of servers. 



in order to satisfy the m-guillotine property, witliout increasing tlie 
total edge length too much. 

We then apply a dynamic programming algorithm, running in 
0{n'^^"'') time, to compute a minimum-cost covering network hav- 
ing a prescribed set of properties: (1) it satisfies the ra-guillotine 
property (with respect to 6), which is necessary for the dynamic 
program to have the claimed efficiency; (2) its disks cover the clients 
Y; and (3) its edge set contains an Eulerian subgraph. This third 
condition allows us to extract a tour in the end. In the proof of the 
following theorem (see the full paper), we give the details of the 
dynamic programming algorithm that yields: 

Theorem 19. The min-cost covering tour problem has a PTAS 
that runs in time 0(n'''''^'). 
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